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Abstract: We study non-Gaussianity in two distinct models of preheating: instant 
and tachyonic. In instant preheating non-Gaussianity is sourced by the local terms 
generated through the coupled perturbations of the two scalar fields. We find that 
the non-Gaussianity parameter is given by f^ L ~ 2g < 0(1), where g is a coupling 
constant, so that instant preheating is unlikely to be constrained by WMAP or 
Planck. In the case of tachyonic preheating non-Gaussianity arises solely from the 
instability of the tachyon matter and is found to be large. We find that for single 
field inflation the present WMAP data implies a bound V^ 4 /Mp < 10~ 4 on the scale 
of tachyonic instability. We argue that the tachyonic preheating limits are useful also 
for string-motivated inflationary models. 
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1. Introduction 

Preheating, first realized in and worked out in detail in Refs. [§, ||, |J, is an 
interesting possibility for obtaining a thermal Universe after the end of inflation. 
Preheating is a non-perturbative phenomenon and perhaps the most efficient way 
of transferring the inflaton energy density into other degrees of freedom. Without 
preheating the inflaton would have to decay perturbatively and slowly, thereby re- 
quiring a long time scale for the thermalization of the inflaton decay products (for 
recent discussions, see 0). 

The simplest realization of preheating is obtained if the inflaton condensate has a 
coupling g 2 (p 2 a 2 , where ip is the inflaton and a is another scalar field. Then during the 
coherent oscillations of ip a resonant production of a quanta can take place due to a 
temporary vacuum instability. The occupation number of a field increases gradually 
while the fluctuations in a increase exponentially. An important observation is that 
during this phase gravitational (metric) fluctuations also get amplified at super- 
Hubble scales 0, [7] (for a general treatment of first order linear perturbation theory, 



The growth of metric fluctuations is of great interest because of the potential 
observational consequences. In particular, during inflation the growth of the second 
order metric fluctuations leads to an enhancement of the non-Gaussianity parameter 
fff L of the CMB (for a review of non-Gaussianity, see In our recent paper [|T^] we 



see §) 
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considered a particular preheating scenario where in comparison with the inflaton, 
the dynamics of a was assumed to have negligible impact on metric fluctuations. We 
estimated the second order metric fluctuations following [11], |T2j and found that ff iL 
grows exponentially with a rate that depends on the number of inflaton oscillations. 
(This is a generic result and applies also to the case of supersymmetric flat directions, 
which could provide an alternative scenario for reheating the Universe [0). The 
result suggests that non-Gaussianity may provide an important test of preheating 
scenario, and in the present paper we shall demonstrate that this is indeed so. 



Non-Gaussianities can also arise after inflation Jm) when the energy in the in- 
flaton condensate is transferred to other degrees of freedom. Such a situation arises 
during preheating as we shall now argue. Let us focus on two different cases. First, 
let us note that the inflaton could decay non-perturbatively during just one oscilla- 
tion only. This is simple to understand for example in a chaotic type model with 
V ~ m 2 (f 2 /2. There inflation ends at <p ~ Mp and the field rolls down towards its 
global minimum. However, by virtue of the coupling, g 2 ip 2 o~ 2 , the a field obtains an 
effective mass m 2 = g 2 <p 2 . Initially, as long as |m CT | <C m 2 , the fluctuations in a evolve 
smoothly. Eventually the adiabatic condition is violated, and particle production oc- 
curs when \m a \ ~ g<p > g 2 tp 2 . During the first oscillation the Hubble expansion rate 
can be neglected, i.e. > H(t). Hence the field velocity near the minimum of 
the potential is given by \<p\ ~ m^ipe, where (p e ~ O.IMp is the amplitude of the 
first oscillation after the end of inflation. In this situation non-adiabatic condition 
is violated when — <p* < (p < <p*, where <^>* ~ \frn^p^j g |3|, [L5|]. For a sufficiently 
large coupling, g ^> 10~ 4 , which is required for an efficient reheating, the interval of 
non-adiabatic violation is very short. This situation is known as instant preheating 



and has many interesting cosmo logical consequences (for a discussion, see [15]). 

The second case is obtained in models where the scalar field undergoes insta- 
bility due to the appearance of a tachyonic mass term — m 2 x 2 /2 in the effective 
potential [|I^| which thus induces symmetry breaking, where we assume x is an_ 
other scalar field besides inflaton ip. In our case we may assume that the symmetry 
breaking occurs smoothly and the tachyonic mass lasts for a very short period much 
smaller than the phase transition time scale and the Hubble rate. During this period 
the long wavelength fluctuations with momentum k < m grow exponentially with 



$Xk ~ exp(V^ 2 — k 2 t). This situation is known as tachyonic preheating. 

Our main aim is to estimate the second order metric fluctuations during these 
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short periods of particle creation occurring in the two preheating mechanisms. En- 
hancement of non-Gaussianity is expected to take place because of the presence of 
second order terms sourced by the first order matter perturbations. For instant 
and tachyonic preheating we need to adopt the multi-field formalism developed in 
Refs. Jl2|, in order to account for non-Gaussianity due to the growth of pertur- 
bations in a and %■ However as we shall see in that tachyonic case we can simplify 
the analysis by assuming that the inflaton VEV vanishes after inflation, which allows 
the metric fluctuations to be mainly seeded by the tachyonic instability. 



2. Basic equations 

Let us first recapitulate the basic equations required for both instant and tachyonic 



preheating. Following Ref. |TTJ the metric decomposition is given by, 

Soo = -a( V ) 2 (1 + 20« + 0( 2 >) , (2.1) 
0m = O, (2.2) 
g ij = a(r J ) 2 (l-2V (1) -V (2) Hi> (2-3) 

where the generalized longitudinal gauge is used and the vector and tensor pertur- 
bations are neglected. Here rj denotes the conformal time and a(rj) the scale factor. 

In the case of instant preheating we need to study the background equations of 
motion for the two fields ip, a. For simplicity and for the sake of clarity we assume 
that the VEV of a vanishes, (a) = 0, which makes it possible to obtain analytic 
approximations from the second order perturbation equations, as we have shown 
earlier [H|. Such a situation occurs if a is driven to the minimum of its potential 
right after the end of inflation due to a positive Hubble-induced mass correction, 
which arises very naturally in many supersymmetric theories (for a review, see ||17|| ). 
For our purposes it is the background motion of (p rather than that of a that is 
important in the instant preheating case. However, it is interesting to investigate 
the fluctuations 5a generated via the coupling to (p. 

In the tachyonic case for simplicity we will assume an opposite scenario, where 
after inflation the inflaton VEV vanishes, (<p) = 0, while the tachyon field x is rolling 
down the potential. This will allows us to estimate the second order metric perturba- 
tions |10[ . Therefore while describing tachyonic instability the same equations apply 
if we replace in the inflaton with the tachyon in the inflaton equations, ip — > Xi an d 
the second field with the inflaton in the equations for the second field, a — > (p. 
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The fields can be divided into the background and perturbation: 

(p = Mv) + <$ (1 V(?7, x) + i<f< 2 Vfa, x) , (2.4) 
a^S^a( V ,^) + Vcr(r/,x), (2.5) 
where the background value for a is assumed to vanish. The background equations 



of motion during preheating are then found [[IT], [12|] to be 

3H 2 = y^ 2 + ~K 2 aV(^) , (2.6) 
= <p'i + m<p' b + a*V'(<p b ), (2.7) 

while the avequation is trivial. Here 7i denotes the Hubble expansion rate expressed 
in conformal time. 



The relevant first order perturbation equations can be written in the form [12 



0(D " _ ^0(1) + 2 (n-^j 4F> ' + 2 f W - ^H^j 0« = , (2.8) 

<5 ( V + 2W ( V - did l 5 {1) a + g 2 <p 2 = . (2.9) 

All the information regarding S^ip is contained in Eq. ( |2.8| ), whose right hand side 
is zero by virtue of (a) = 0. Further note that there are no metric perturbations in 
Eq. Q2.9p . This is due to assuming a vanishing VEV for a. Now the a part can be 
solved separately and for the rest the usual single-field results [jTTJ apply. 

At the second order we are only interested in the gravitational perturbation, 



whose equation can be written in an expanding background as [12 



+2^n-^j ^' +2(h'- ^-nj ^ - ^v (2) = 

local local i^non-local j 

(2.10) 

where the source terms J are quadratic combinations of first order perturbations 



<Jtp, local rt 



-2(^V) 2 - 8(^) 2 (0 (1) ) 2 + 8^</WV + a 2 f^ (1 V) 



d(p 2 



-24H'(f') 2 - 24H<i> w <i> w (2.11) 
X,iocai = -2n 2 (6^ar + M^(6®<t? , (2.12) 

Jnon-local = A" 1 /^^, ^ ^ ^ , (2-13) 

where / is a quadratic function of the first order fluctuations and the coefficients de- 
pend on background quantities. Because of the inverse Laplacian the last source term 
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is non-local. Typically such term contains: A A(f)^), A 1 (di5^cpd i 5^(p), . . . 

Note that the left hand side of Eq. (|2.10|) is identical to the first order equation, see 
Eq. (£§)■ 



v = ^ m y + gV*\ (3.1) 



3. f$[ L for Instant Preheating 

Let us now consider a two-field model with the potential 

1 

—1 
2 

where (p is the inflaton condensate with a mass m v . In instant preheating the particle 
production occurs during one oscillation of the inflaton. The particle production 
occurs when the inflaton passes through the minimum of the potential ip = 0. In 
this case the process can be approximated by writing 

up = <p (t - t ) , (3.2) 

where ipo is the velocity of the field when it passes through the minimum of the 
potential at time to. The time interval within which the production of a quanta 
occurs is 

At« = (<^o|)- 1/2 , (3-3) 

which is much smaller than the Hubble expansion rate; thus expansion can be ne- 
glected. Note that by virtue of the coupling g 2 ip 2 a 2 the a field acquires a mass and 
provided that g < H inf /(p ~ 1(T 5 for H inf ~ 10 13 GeV and <p ~ k' 1 ~ 10 18 GeV, 
the fluctuations in a field were already present on large scales during inflation with 
8®o ~ H mf /2n. 

The occupation number of produced particles jumps from its initial value zero 
to a non-zero value during — ip* < ip < tp*. In the momentum space the occupation 



number is given by [15 



/ \ 

n k = exp — - , (3.4) 

V 9\<Po\J 

and the largest number density of produced particles in x-space reads 

(3.5) 

with the particles having a typical energy of (g|0 o |/7i") 1//2 , so ^ na ^ their total energy 
density is given by 

l a (i)-N 2 (ffl^ol) 2 f o R \ 

Pa~^ 6 V) ~ g7r 7/2 • ( 3 ' 6 ) 
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These expressions are valid if m 2 < g\<fio\, a condition that we assume for the rest of 
our calculation. 

Ignoring the expansion of the Universe (a = 1, rj = t), and using Eq. (|3.2|), the 
second order gravitational perturbation, Eq. (|2.1Up , is at large scales 

$ 2 >~-2« 2 (<f«<7) 2 . (3.7) 

The non-Gaussianity in the gravitational potential, parameterized as = /^(0 ( ' 1 ^) 2 , 
can be estimated by solving the second order gravitational potential from Eq. ( |8.7| ) 
using Eq. ( |3.6|) . We obtain 



NL 



^ 2 |0o| 2 At 2 



(3.8) 



?7T 7 /2Af|,(0(l))2 

It is a simple exercise to estimate the right hand side for a chaotic inflaton potential 
with m = 10 13 GeV and (p^ ~ 10~ 5 . The velocity of the scalar field at the potential 
minimum comes out to be \(fio\ ~ 10~ 7 Mj,; using these values and Eq. ( J3.3|) we 
obtain an estimate for the upper limit of the non-Gaussianity parameter in the case 
of instant preheating: 

fL~29- (3-9) 

Efficient reheating requires g > 10~ 4 . However, since Eq. (|3.9|) implies that ffj L <C 
0(1), instant preheating is unlikely to yield any detectable non-Gaussian signal in 
the forthcoming CMB experiments. The lowest observable value for ffj L by WMAP, 
Planck and an ideal experiment is respectively 13.3; 4.7; and 3.5; including polariza- 



tion data decreases the limits respectively to 10.9; 2.9; and 1.6 [18 



4. ffc L for Tachyonic Preheating 

In order to understand the non-Gaussianity triggered by the tachyonic instability, let 
us assume a simple toy model where there is an inflationary sector and a symmetry 
breaking phase transition with a mass squared term as negative: 

V = V( ( p) + Vo-^m 2 x 2 + j X 4 - (4-1) 

We assume that the inflaton potential is some polynomial potential with a vanishing 
VEV, V(cp) ~ /(y? n ). Inflation is supported by both V((p) + Vo- During inflation 
we assume that the tachyon field is sitting at the maximum x = by virtue of 
large friction. The mass of x is such that the tachyonic instability is triggered when 
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m > H ~ Vo/3Mp. During this period we assume that the inflaton settles down 
to ((f) = 0. This will allow us to separate the tachyon fluctuations from that of 
the inflaton. This also allows us to use the same equations ( |2.4| - 12.13| ) but now the 
tachyon field x obeys Eq. fl2.8|) and the inflaton field ip obeys Eq. (^.9|). So we replace 
ip — > x i n Eq. (|2.8| ) and in a — > if in Eq. ( |2.9|) , and we make similar replacements in 
Eq. (|2.10|) and the expressions for J x .\ oc& \ and ^locai respectively. 

The rolling of a tachyon in itself results in an exponential instability in the 
perturbations of x with physical momenta smaller than the mass. The tachyonic 
growth takes place within a short time interval, £* ~ (l/2m) ln(7r 2 /A) (see ||16|| ). 
During this short period the occupation number of x quanta grows exponentially for 
modes k < m up to ~ exp^mt*) ~ exp(ln(7r 2 /A)) ~ 7r 2 /A. For very small self- 
coupling, which is required for a successful inflation, the occupation number, which 
depends inversely on the coupling constant, can become much larger than one. 

The scalar field fluctuations, which are responsible for exponentially enhancing 
the occupation number for x quanta, also couple to the metric fluctuations. If we 
assume that the modes grow within a time interval much smaller than the Hubble 
rate, we can set TC = in Eq. ( |2.8| ). Then, in the long wavelength limit, we get from 
Eq. @, 

-2A0 (1) = 0, (4.2) 

where A = Xo/Xo- With the assumption of a brief tachyonic stage, x , Xo are ef- 
fectively constants. Note that although during rolling tachyon the long wavelength 
modes are excited, but it is important that the tachyon perturbations must exist 
during inflation. In this respect the tachyon fluctuations are isocurvature in nature. 
In order to further simplify our calculation we neglect the inflaton perturbations in 
our subsequent analysis. 

There are two solutions of Eq. ( f4.2| ); a constant (p^ ~ 10~ 5 , and an exponentially 
growing solution <fy$ cx exp(2 At). The former case arises when we recognize the 
constant by the temperature anisotropy of the observed CMB fluctuations. If the 
isocurvature component at the end of inflation is small, then the first derivative of 
</>W is also small but non-vanishing. Hence we may neglect the exponential solution 
of the first order metric perturbation. Although our argument holds good for the 
first order metric perturbations, but as we shall show this will not be the case for the 
second order calculation. With these simplified approximations we can then estimate 
the amount of generated non-Gaussianity by following a logic similar to the case of 
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instant preheating. 

First, the number density of the produced particles in x-space is given by n x ~ 
m 3 /(87rA). Hence the total energy density stored in produced x quanta is given by 

p x ~i(5«x) 2 ~mn x ~^. (4.3) 

The equation for the second order metric perturbation now includes a source term 
which includes i7r,iocai an d .TxJocal- However the main contribution comes from the 
excitations of the tachyonic instability. The inflaton fluctuations are subdominant 
compared to the exponential growth of when it is settled around its VEV 

((p) = 0. In the the long wavelength regime the perturbation equation reads as 

fiH^-l^. (44) 
7T A 

Integrating the above equation over the time interval ~ (l/2m) ln(7r 2 /A), we find 
0( 2 ) r*j {m/Mp) 2 ln 2 (7r 2 /A)/(47rA). The non-Gaussianity parameter for tachyonic pre- 
heating in case the first order metric perturbation stays constant is then roughly 
given by 

St,^U^~M(^), (4-5) 



4n\M P J A(0«) 2 V A 
where we substitute k ~ Mp l . Writing this in terms of Vq = m 4 /(4A) and taking 
Si 1 ) ~ 10~ 5 , we obtain 



f* L ~ 1.6 x 10 9 A-/ 2 i^pj In 2 . (4.6) 

This expression should be compared with the observationally constrained one: f^L — 
~ftiL + see for instance H]. 1 The current observational limit of the non- 



linearity parameter set by WMAP is —132 < f^ L < 60, at 95% confidence level [[IS 

Adopting the upper limit \fff L \ < 132 and rearranging, we arrive at the bound 

V^/Mp < 3 x lO^A 1 / 4 In" 1 (tt 2 /A) . 

For an effective field theory to remain perturbative we should require that A <C 1, 

which yields the interesting constraint Vq /Mp <C 10 -4 . Note that compared to the 

1 The measure of non-Gaussianity is the non-linearity parameter (f>i 2 ) = ftiL{4>^) 2 ■ In general 
ffj L contains momentum dependent part, i.e. /^ L (ki,k2), and the constant piece. It is the non- 
local terms which affect the momentum dependent part, since all the derivatives are replaced by 
momenta in the Fourier space. However the present constraint on non-Gaussianity parameter from 
WMAP does not give the momentum dependent constraint but only the constant part. Therefore 
the non-local terms do not lead to any observable constraints, so we do not consider them here. 
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log 10 (A.) 



Figure 1: The parameter space in tachyonic preheating allowed by the WMAP data 
for single field inflation (below the red curve). Black horizontal line is the usual limit 
= 10 16 GeV coming from the COBE normalization. 

usual bound Vq^ < 10 16 GeV from COBE normalization, the absence of observable 
non-Gaussianity implies a bound on the scale of tachyonic instability V^ A which is 
more stringent by two orders of magnitude. The parameter space allowed by WMAP 
data is given by the region below the red curve in Figure [I]. 

When obtaining the result above we have assumed that the first order metric 
perturbations is roughly given by the constant value determined by the inflationary 
epoch. Let us investigate the other limit when the first order metric fluctuations 
also obtain an exponentially growing solution by virtue of the tachyon excitations. 
To check this possibility, let us assume that the exponential solution actually domi- 
nates over the constant one. Following the second order analysis and assuming that 
the main contribution to the second order perturbation arises from the tachyonic 
instability, we obtain from Eq. ( |2.10|) 



4>W _ 2A ^) = _«2 [ 2 (5«x) 2 + 8^(0 (1) ) 2 - V XX (6^ X ) 2 ~ 8xo0 (1) 5 (1) x] , (4.7) 



where V 



m 2 + 3 A% 2 and where can be solved through the Einstein 
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For the tachyonic region 3Ax 2 < rn 2 so that we can take V xx ~ — m 2 . Now 
< - 2 - ) contains the homogeneous solution ~ exp(2At) together with a source part 
' exp(4 At) 2 . After a while the source part dominates and we obtain the result 3 



NL 



(1)^2 



2m 2 4x 2 ^Xo 



K 2 xl k2 Xo xl 



(4- 



The tachyonic growth persists until x ~ m l (2\/A), for a time span £* ~ (l/2m) ln(7r 2 / A) 
[|Tj|. With these approximations, writing V = m 4 /(4A), we obtain 



XMp 1 / V XM P , 



4 /7r 2 



If we assume COBE normalization V^ 4 < 10 16 GeV, the minimum of Vq is given 



by the conditions V Q 1/4 /M P = A 1 / 4 ^ - f$ L and A = 7r 2 exp(-^/(8 - /^ L ) 2 /2 - 8). 
These imply the limit ff iL < —37 regardless of the value of A, well within the 
observational capabilities of WMAP. 

Many string-motivated inflationary models could thus be constrained by the 
present and future limits on non-Gaussianity. An example would be inflation first 
driven by brane-anti-brane interaction and then coming to an end when the tachyonic 



instability is triggered [20]. Should we take the tachyon coupling to be very small 
A ~ 10~ 12 , as constrained by the amplitude of the CMB scalar fluctuations ]2"I], we 
immediately obtain from Eq. (|4.6| ) that V^ 4 /Mp < 10~ 7 in order to comply with the 
current WMAP limit on non-Gaussianity. Such considerations imply very interesting 
constraints on the scale of the tachyonic instability and on the tachyon self coupling 
in brane-anti-brane driven inflation. 

In obtaining our limits we have made some approximations. In particular, we 
assumed that the VEV of a field is vanishing in the instant preheating case and in 
tachyonic preheating we assumed that the inflaton VEV is vanishing after the end of 
inflation. These assumptions make our analysis simple and provide a handle on the 
non-Gaussianity parameter. 



2 The second order metric perturbations always have a growing source term by virtue of the 
non-vanishing background motion of the scalar field, i.e. S^x, S^X- 

3 Assuming A is constant; in reality there will be a small time variation but we may assume that 
most of the interesting modes are growing within a time interval which is short compared to the 
variation in Xo/xo and the Hubble rate. 
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